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INTRODUCTION 


Accurate  drag  predictions  are  essential  for  optimum  missile 
. design.  For  Mach  numbers  in  the  transonic  regime,  the  NSWC  code^ 

is  currently  applicable  to  missiles  with  pointed  or  slightly 

blunted  noses.  This  restriction  results  from  the  use  of  methods 

' 2 

‘ based  on  small  perturbation  theory  such  as  Wu  and  Aoyoma.  When 

f bluntness  effects  become  large,  it  is  necessary  to  account  for 

! radial  flow  velocities  comparable  to  free  stream  velocities  over 

a sizable  region  of  the  flow  field.  A fundamental  failing  of 
small-disturbance  theory  is  that  it  cannot  do  this. 

A precise  calculation  of  nose  wave  drag  depends  upon  the  accu- 
rate determination  of  surface  pressures  both  in  the  stagnation  region 
and  the  shoulder  region.  This  is  because  the  overpressure  in  the 
stagnation  region  creates  a large  drag  force  which  is  partially  or 
wholly  counterbalanced  by  the  flow  expansion  in  the  shoulder  region. 
In  subsonic  flows,  the  balance  is  exact  for  the  nose  itself,  while 
in  transonic  flows,  one  is  attempting  to  find  the  difference  of  two 
large,  nearly  equal  quantities  in  order  to  compute  the  nose  wave 
drag.  Consequently,  a method  is  required  which  can  maintain  the 
required  accuracy  in  these  regions. 

The  computation  of  nose  wave  drag,  then,  for  the  transonic 

, speed  range  requires  an  accurate  code  capable  of  properly  treating 

‘ the  inviscid  flow  physics.  The  exact  governing  equations  are  the 

'•  Euler  equations.  In  the  present  report  an  implicit  Euler  equation 

* 3 4 5 

solver  ' ' is  used  to  calculate  the  flow  field  around  the  bodies 

y 

f and  pressures  on  the  bodies.  These  pressures  are  then  integrated 

( to  obtain  the  nose  drag  coefficient  for  a set  of  Mach  numbers,  nose 

1 bluntnesses  and  nose  lengths.  The  resulting  table  of  transonic 

[ drag  coefficients  which  is  necessarily  restricted  to  sphere-tangent 

: ogive-noses  is  incorporated  into  the  NSWC  computer  code.^ 
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ANALYSIS 


Implicit  Unsteady  Euler  Equation  Method 

The  partial  differential  equations  governing  the  two-dimensional 
planar  or  axisymmetric  flow  of  an  unsteady  inviscid,  non-heat-con- 
ducting,  ideal  gas  are  written  in  conservation-law  form  under  the 
generalized  independent  variable  transformation 

T = t, 

C = C (t,x, y) , 
n = n(t,x,y) 

as  follows: 

(U/J)^  + [(C^U  + CjjE  + C^F)/J]^  + (n^u  + + nyF)/j]^  + jH  = 0 (1) 

where 


pv 
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PU 

2 
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and  the  Jacobian 


J 


^x'^y 


In  eq.  (1) , j = 0 for  two-dimensional  flow;  j = 1 for  axisym- 
metric flow  (where  y becomes  the  cylindrical  radius) ; p represents 
the  pressure,  p,  the  density;  u and  v,  the  velocity  components  in 
the  X and  y directions,  respectively;  and  e,  the  total  energy  per 
unit  volume.  The  following  equation  relates  the  pressure,  density, 
and  velocity  components  to  the  energy  for  an  ideal  gas 

e = p/(Y-l)  + p(u^  + v^)/2  (2) 


The  equations  represented  by  eq.  (1)  were  solved  using  an  implicit 
Euler  equation  m.ethod. 
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Using  the  above  independent  variable  transformation,  the 
physical  plane  with  boundaries  1-4  in  x,  y,  and  t space  is  mapped 
into  the  rectangular  computational  plane  n,  and  t space  as 
shown  in  figure  1.  The  metrics  (geometric  factors)  that  are 
required  by  eq.  (1)  are  defined  as  follows: 


’t ' <y{*T  - 


''x  ■ 

n = x./i 

y C 


(3) 


where 


I = l/J  = Xjy^  - x^y^. 

Generally  at  each  time  step  the  metrics  of  eq.  (3)  are  not  known 
analytically  and  they  must  be  determined  numerically. 

The  basic  numerical  algorithm  used  in  the  IMP  code  was  developed 

4 

by  Beam  and  Warming.  It  is  second-order  accurate,  noniterative, 
and  is  in  a spatially  factored  form.  In  particular,  the  "delta-form" 
with  Euler  time  differencing  is  used,  and  when  applied  to  eq.  (1) , 
it  assumes  the  following  form 


.n. 


.n . 


,n , 


(4) 


(I+Atd^A’^)  (I+At6  b")  = -At(6^E“+6  F“+Hj) 

C n C n 

where  A and  B are  the  Jacobian  matrices  A = 8E/9u  and  B = 9F/9u 
with  u = U/J;  E = + 5yF)/J,  and  F = (n^U  + + nyF)/J. 

The  A matrix  is  given  below.  The  B matrix,  which  has  the  same  form 
except  that  ^ ^ , and  C are  replaced  by  nx.f  n f and  n , is  not 

t X y u X y 

presented.  I is  the  identity  matrix,  6c.  and  6 are  second-order, 

/vn+l  z'  ^ 

central  difference  operators,  U = U(nAt)  and  At  is  the  integration 
stepsize. 
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rj  = 77(1, x,y) 


Figure  1.  Transformation  of  Physical  Plane  to  Computational 

Plane  for  a Single  Solver 
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Equation  (4)  is  solved  at  the  interior  points  only.  It 
requires  two  block  4x4,  tridiagonal  inversions  at  each  time  step 
of  the  integration.  The  solution  procedure  is  as  follows; 

A n4'  1 ^ n 

1.  Define  AU  = U - U ; 


2. 


3. 


4. 


5. 


6. 


Form  the  right-hand  side  of  (4)  and  store  results  in  the 
array; 


'^n+1  *^11+1 

Apply  smoothing  U = U - e/8  S/J  where  the  smoothing 
term  is  defined  below; 


Define  U = (I  + At6pB*')  U and  solve  the  matrix  equation 
(I  + At6^A”)U  = for  U and  again  store  the  result  in 

the  array; 


Tl  ^ ^ 

Solve  the  matrix  equation  (I  + At6^B  ) AU  = U for  AU; 
Obtain  the  new  values  of  U from  the  equation 

An+l  ^ 

V ^ = AU  + U . 


i 


M 


5 


A fourth  order  smoothing  term  is  needed  to  eliminate  nonlinear 
instabilities^  which  arise  since  the  use  of  central  differences  in 
the  spatial  directions  results  in  a neutrally  stable  algorithm. 

J 

The  smoothing  term  S is  | 


I 

1 

J 


I and  C is  the  smoothing  coefficient  chosen  from  the  range  OieiO.4 

i depending  upon  the  size  of  the  time  step.  The  j and  k indices 

j correspond  to  the  C and  n directions,  respectively.  At  the  points 

adjacent  to  the  boundaries  a special  form  of  the  smoothing  term  is 

I 

used. 

The  solution  of  initial  value  problems  using  finite  difference 
procedures  requires  an  integration  step  size.  The  value  of  this 
step  size  is  a function  of  the  eigenvalues  of  the  coefficient 
matrices,  A and  B,  of  the  governing  pde's  (equation  1)  when  written 
in  the  form; 

u + AUp  + Bu  + jH  = 0 (7a) 

T 5 n 


n+1 


n+1 


n+1 


, n+1 


n+1 


,A,,n+l 


The  eigenvalues  of  the  coefficient  matrices  A and  B are 


K = C,.  + + vC  ± c(5^  + 

=54-  + + v£ 


1/2 


’b,  = ^t  ^ 

X ^ A 


1/2 


'B, 


= Tit  + ^^x  + 


where  H is  the  modified  source  term  and 


(7b) 


1/2 

where  c,  the  speed  of  sound  is  equal  to  (Yp/p)  . The  a terms  in 
Eq.  (7b)  can  be  recognized  as  the  slopes  of  the  characteristics  and 
streamlines  in  the  t-C  and  t-n  planes. 

The  stability  theory  for  explicit  methods,  at  least  for  (t,C) 
space  in  unsteady  flow,  requires  that 

f|  i . 1,3  (8b) 

' 1(0*)  I 

Z max 

where  (o,  ) is  the  maxima  of  the  local  eigenvalues  of  all  the 
max 

points  at  a particular  constant  t or  5 plane.  Similar  relations 
can  be  obtained  for  the  (t,n)  plane: 

At  X 1 

An  Z = 1,3  (8c) 

(Ob  ) 

Z max 

This  planar  analysis  has  been  shown  to  give  a good  bound  on  the  step 
size  in  multidimensional  problems  if  the  right  hand  side  of  Eqs. 

(8b)  and  (8c)  are  multiplied  by  a constant  C.  In  explicit  methods, 

C 1 and  is  usually  assigned  a value  of  0.9  whereas  in  implicit 
methods  C is  only  restricted  by  the  accuracy  of  the  solution  since 
implicit  methods  are  unconditionally  stable.  Typically  these  values 
are  on  the  order  of  5 to  10  for  inviscid  flow  computations. 

It  is  convenient  to  consider  two  types  of  boundaries  in  inviscid 
computational  fluid  dynamics  - those  which  are  impermeable  (no  mass 
flow  across)  and  those  which  are  permeable  (mass  flow  across) . The 
impermeable  boundaries  consist  of  such  surfaces  as  solid  walls, 
symmetry  planes  and  slip  surfaces,  while  the  permeable  boundaries 
consist  of  such  things  as  supersonic  inflow  or  outflow  boundaries 
and  shock  waves.  In  the  transonic  flow  about  projectiles,  the  four 
boundaries  of  figure  1 can  be  characterized  as  follows:  boundary  1 
is  either  a shock  wave  (M^  > 1.00)  or  a fixed  outer  boundary  (M^<1.0), 
boundary  2 is  an  outflow  plane,  boundary  3 is  the  solid  body,  and 
boundary  4 the  plane  of  symmetry. 
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The  body  boundary  is  treated  by  using  the  tangency  condition 
which  is  incorporated  through  the  use  of  Kentzer's  scheme,®  and  at 
the  symmetry  plane,  the  variables  are  reflected  according  to  whether 
they  are  odd  or  even.  The  outflow  boundary  is  placed  far  enough 
down  the  aft  body  such  that  the  pressure  can  be  approximated  by  its 
value  at  infinity  and  the  rest  of  the  variables  are  extrapolated. 

At  the  outer  boundary,  there  are  two  possibilities:  (1)  if  shock 
fitting  is  used  and  the  flow  is  supersonic,  the  sharp  discontinuity 
approach  of  Thomas^  is  used  or  (2)  if  shock  capturing  is  used,  the 
flow  can  be  either  subsonic  or  low  supersonic  and  the  outer  boundary 
is  chosen  such  that  its  influence  on  the  flow  field  surrounding  the 
body  is  negligible. 

The  initial  conditions  are  determined  by  first  guessing  a shock 
wave  shape  or  outer  boundary  shape  for  the  particular  obstacle  at 
hand.  Second,  a Newtonian  pressure  distribution  is  prescribed  on  the 
body.  Using  the  fact  that  there  is  no  flow  through  the  body  and  that 
the  maximum  entropy  streamline  wets  the  body,  the  remainder  of  the 
properties  can  be  prescribed  on  the  body.  Finally,  a linear  variation 
for  the  properties  is  prescribed  between  boundary  1 and  the  body. 

This  gives  an  initial  guess  for  the  flow  field  which  is  integrated 
forward  in  time  until  a steady  state  is  reached. 

In  all  cases,  the  resulting  flow  field  between  the  body  and  the 
outermost  boundary  is  treated  in  shock  capturing  fashion  and,  therefore, 
allows  for  the  correct  formation  of  secondary  internal  shocks. 

Applying  the  algorithm,  boundary  conditions  and  initial  con- 
ditions to  Eq.  (1)  and  solving  Eg.  (1)  in  a time  asymptotic  fashion 
with  interest  only  in  the  steady  state,  a complete  set  of  flow  field 
data  results  out  of  which  the  nose  drag  is  calculated  by  integrating 
the  surface  pressures  over  the  surface  area  of  the  body. 

Drag  Calculation 

The  wave  or  pressure  drag  for  a sphere-ogive-cylinder  body  at 
zero  angle  of  attack  is 
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D = J J (P-P„)y  dy  d(}) 

0 0 

(see  figure  2) . Note  that  p-p„  is  the  differential  pressure  along 
the  surface  and  yd<{>*dy  is  the  surface  area  over  which  the  pressure 
is  acting  for  drag  considerations.  Since  the  body  is  axisynunetric 
and  at  zero  angle  of  attack,  the  pressure  is  not  a function  of  the 
meridional  angle,  (f).  Integrating  the  above  equation  (9)  and  using 
the  fact  that  p-p„  = f(y)  gives 

^b 

D = 2-n  (p-p<„)ydy  . 
j 
0 


Figure  2.  Nomenclature  for  Drag  Calculation 
By  using  the  following  approximation  for  total  drag 


equation  (10)  can  be  integrated  using  the  mean  value  theorem  for 
each  geometric  segment  i.  This  can  be  written  as  a form  of  the 
trapezoidal  rule 

2 2 

I (y:  - yT.j) 

D 2Tr  I (p-p^)  2 ^ 

i=l  ^ 


(10c) 


where 


(P-Po,)  = 


(P-P„)i  + (P-Poo)i_i 


Thus  the  total  drag  is  the  summation  of  each  geometric  increment  of 
drag. 

The  drag  coefficient  is  defined  as 


where 

1 2 

q is  the  dynamic  pressure  -jPooV^ 

is  the  reference  area  irrj^ 

Finally,  the  total  drag  coefficient  is 


iniax  IS-  ^ Sj.J  2 „2  . 

2 . ^ 2 ^i  ^i-1 

^b  ^ 


(lib) 


where 


P-P„ 
C = 

p q 


Body  Geometry 

A large  number  of  ogive  body  shapes  are  considered  in  this 
report.  They  are  characterized  by  bluntness  ratios  from  0 to  1 
and  fineness  ratios  from  5 to  0.5. 

From  figure  3,  the  bluntness  ratio  is  defined  as  the  nose  radius 
divided  by  the  cylindrical  body  radius 

Rn 

Bluntness  = ^ 

% 


i 


(0,0)[  "" 


IJ. 


Figure  3.  Ogive  body  geometry  definition 


where  a bluntness  of  0 is  a pointed  ogive  and  a bluntness  of  1 is  a 
hemisphere.  The  fineness  ratio  is  defined  as  the  sphere  ogive  body 
length  divided  by  the  diameter  of  the  attached  cylinder: 

L 

Fineness  = 

where  a fineness  ratio  of  0.5  is  a hemisphere  cylinder. 

Geometrically,  the  sphere  ogive  cylinder  configuration  is 
defined  by  using  one  algebraic  equation  in  y and  x to  define  each 
of  the  three  segments  which  are  shown  in  figure  3.  This  general 
equation  is  written 

2 2 

Ay  + Byx  + Cx  + Dy  + Ex  + F = 0 

where  the  coefficients  vary  depending  upon  which  segment  of  the  body 

is  being  defined.  Body  segment  1 is  input  as  the  arc  of  a circle 

of  radius  with  the  origin  at  and  it  satisfies  the  equation 

y^  + (x-R  )^  = R^  . 

n n 
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similarly,  body  segment  2 is  also  an  arc  of  a circle  of  radius  R2 
with  the  origin  at  [L^,  (Rj^-R^)],  and  satisfies  the  equation 


The  third  segment  is  a straight  line  parallel  to  the  x axis  and  is 
given  by 

y = Rb  • 

t 

I In  the  computer  solution,  all  of  the  geometric  quantities  are 

I 

j normalized  with  respect  to  the  cylindrical  radius  Rj^. 

■ Table  Lookup  Method 

1 

j The  wave  or  pressure  drag  coefficient  has  been  calculated  as 

a function  of  three  variables, 

Cp  = f(M.,  R„/R^,  VRj,)  . 

The  problem  is  one  of  interpolation  for  the  drag  coefficient  given 
the  three  independent  variables.  This  can  be  done  by  three  one- 
dimensional interpolations.  To  keep  within  the  framework  of  the 
previous  program  usage  and  for  the  ease  of  understanding,  the 
following  approach  was  taken.  In  this  process,  the  table  of  drag 
coefficients  consists  of  10  nose  bluntnesses,  8 nose  lengths,  and 

! 4 Mach  numbers.  The  procedure  is  as  follows: 

1.  For  a single  value  of  bluntness,  load  up  a new  set  of  8 
nose  length  arrays  each  consisting  of  drag  values  for 
the  four  corresponding  Mach  numbers;  i.e.  = L^(Cjj(M^)) 

X ~ 1,  2,  ...,  8. 

i 

I 

I 

I 


i 
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2.  For  a given  Mach  number,  a three  point  Lagrangian  inter- 
polation is  used  to  obtain  a value  of  the  drag  coefficient 
from  each  of  the  8 new  nose  length  arrays  of  step  1.  Each 
interpolated  value  is  stored  in  a new  array  corresponding 
to  one  of  the  8 nose  lengths. 

3.  For  a given  nose  length,  a five  point  Lagrangian  interpolation 
is  used  on  the  new  arrays  created  in  step  2 to  find  the 

final  value  of  drag  which  is  stored  in  a new  array. 

4.  Steps  1 through  3 yield  the  drag  coefficient  for  a single 
value  of  bluntness  in  the  drag  table  and  for  the  given 
Mach  number  and  nose  length.  Repeating  steps  1 through  3 
for  each  value  of  nose  bluntness,  a value  for  the  drag 
coefficient  is  obtained  for  each  bluntness. 

5.  A five  point  Lagrangian  interpolation  is  performed  on  the 
drag  coefficient  array  of  step  3 which  consists  of  a 
value  for  each  of  the  10  bluntness  in  the  table.  This 
gives  the  final  value  of  the  drag  coefficient  for  a 
specified  Mach  number,  nose  length,  and  nose  bluntness. 

Two  subroutines  not  available  in  the  NSWC  parent  code  were  needed 
in  this  solution  procedure.  They  are  documented  in  a later  section. 

RESULTS  AND  DISCUSSION 

Comparison  of  Pressure  Prediction  Methods 

0 

Since  there  is  experimental  data  readily  available  through 
the  transonic  range  for  hemisphere  cylinders,  the  present  technique 
was  initially  applied  to  that  limiting  case.  These  results  are 
presented  in  figures  4-10. 

Figure  4 shows  the  typical  grid  network  used.  By  numerical 
experimentation,  it  was  found  that  locating  the  outer  boundary  at 
10  body  diameters  and  the  outflow  boundary  at  4 body  dieuneters  was 
sufficient  to  obtain  a solution  uninfluenced  by  the  choice  of 
boundaries.  It  was  also  found  useful  to  stretch  the  mesh  in  the 
outer  regions  and  compress  it  near  the  body. 
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Figure  5 shows  numerical  results  from  the  present  method 
compared  to  the  experiments  performed  by  Hsieh  for  = .7,  .95, 

1.0,  1.05,  1.1.  In  all  cases,  the  comparison  is  very  good  with 
experiment.  This  establishes  confidence  in  the  present  method. 

Figure  6 compares  numerical  computations  of  the  drag  for  a 
hemisphere-cylinder  with  Hsieh 's  experimental  results.  The  drag 
calculation  by  the  present  method  gives  good  agreement  with  the 

O 

previous  niimerical  and  experimental  results  of  Hsieh. 

The  flow  field  of  a highly  supercritical  flow  is  presented 
in  figure  7.  The  typical  transonic  flow  pattern  is  observed  with 
an  acceleration  of  flow  over  the  hemisphere  nose  to  sonic  velocity 
and  downstream  compression  on  the  cylinder  via  a normal  shock  wave. 

Figures  8 and  9 show,  respectively,  the  shock  angle  which  is 
measured  from  the  horizontal,  6,  and  the  shock  wave  and  sonic  line 
for  low  supersonic  flow,  = 1.1  flow  over  a hemisphere-cylinder. 

In  figure  9,  one  of  the  more  interesting  features  of  the  flow 
field  is  the  bending  back  of  the  sonic  line  past  the  vertical  axis. 

Figure  10  shows  the  drag  (pressure)  evolution  for  transonic 
flows  over  a hemisphere-cylinder.  Note  that  in  general  the  limiting 
value  for  drag  is  reached  by  the  time  1200  iterations  have  been 
performed . 

The  results  presented  above  have  established  the  accuracy  and 
versatility  of  the  present  code.  The  next  step  is  to  apply  it  to 

2 

transonic  flow  over  sphere-ogive-cylinders.  Previous  researchers 
have  used  an  inviscid  transonic  approximation  theory  to  calculate 
the  entire  flow  field  around  nonlifting  slender  ogive  bodies.  The 
transonic  flow  fields  are  essentially  the  same  for  these  shapes 
which  include  the  hemisphere  cylinders  as  a special  case.  Figures 
11  and  12  give  some  results  for  different  sphere-ogive-cylinder 
bodies.  The  calculated  pressure  distribution  along  the  body  for 
= 0.95,  4 caliber  tangent  ogive  body  is  presented  in  figure  11 
along  with  available  experimental  data.  The  final  mesh  point  dis- 
tribution used  in  calculating  the  pressure  distribution  was  used 
in  the  remaining  calculations  based  on  its  ability  to  capture  the 
essentials  of  a flow  field  surrounding  sphere-ogive-cylinders. 
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(c)  M = 1.0 


X/R. 


(d)  = 1.05 

Figure  5.  Concluded 
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Figure  6.  Hemisphere-Cylinder  Drag  Variation  with 
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Figure  7.  Embedded  Supersonic  Region  at  = 0.95  for  Flow 
Over  a Hemisphere  Cylinder 
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Figure  10.  Evolution  of  Drag  Solution  for  He 
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Figure  11.  Comparison  Between  Experiment  and  Theory  for 
M =0.95  Flow  Past  a 4 Caliber  Tangent  Ogive  Body 
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Figure  12.  Calculated  Body  Pressure  on  Sphere-Ogive-Cylinders  for 
Varying  Nose  Bluntnesses,  Nose  Lengths  and  Mach  Numbers 
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Figure  12.  Concluded 


Typical  results  for  transonic  flow  about  sphere-ogive-cylinders 
are  presented  in  figure  12.  Figures  12a  and  12b  show  the  effect  of 
nose  bluntness  and  nose  length  on  the  body  pressure  distribution  for 
a = 1.05.  The  major  change  in  the  low  supersonic  flow  regime 
takes  place  near  the  nose  of  the  vehicle.  In  figure  12c  at  a high 
subsonic  Mach  number,  = 0.95,  there  occurs  two  major  changes, 
one  aft  of  the  ogive-cylinder  juncture  and  one  in  the  region  of  the 
nose  as  in  the  low  supersonic  case.  The  major  change  aft  of  the 
ogive-cylinder  juncture  is  due  to  the  termination  of  the  supersonic 
embedded  region  by  a normal  shock  wave  which  produces  a steepening 
of  the  pressure  gradient.  In  figure  12c,  this  phenomena  occurs  at 
about  1.5  calibers  downstream  of  the  nose.  For  decreasing  nose 
bluntness,  the  location  of  the  normal  shock  is  located  more  aft  on 
the  body.  It  should  be  understood  that  in  figure  12c  the  curve 
for  was  obtained  by  placing  more  points  on  the  afterbody 

where  the  terminator  shock  wave  strikes  the  body.  This  in  no  way 
reflects  upon  the  accuracy  of  the  drag  calculation  since  those  were 
done  by  maximizing  the  number  of  points  in  the  nose  region  where 
tiie  resolution  of  the  pressure  is  important  for  the  calculation  of 
drag. 

Drag  Predictions  and  Comparisons 

The  purpose  of  this  study  was  the  calculation  of  the  nose  drag 
coefficients  for  a family  of  sphere-ogive-noses  through  the  transonic 
Mach  number  range.  A table  was  constructed  (Table  I)  for  these 
coefficients  as  a function  of  nose  bluntness,  nose  length,  and  Mach 
number.  This  table  was  incorporated  into  the  NSWC  computer  code. 

As  a first  step  in  obtaining  the  aforementioned  table,  nose  drag 
coefficients  are  presented  in  figures  13  as  a function  of  nose 
bluntness  and  nose  length  for  the  transonic  Mach  numbers  0.8,  0.95, 
1.05,  and  1.2.  To  reduce  the  overall  computation  time  of  this  study, 
a minimum  matrix  of  variables  was  necessarily  calculated.  These 
computations  are  represented  by  the  symbols  in  figures  13. 

In  figures  13a,  c,  e,  g,  the  nose  drag  coefficient  versus  nose 
length  with  nose  bluntness  as  a variable  is  presented.  For  all 
bluntnesses  and  for  all  Mach  numbers,  the  nose  drag  coefficient 
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TABLE  I.  Cjj  = Sphere-Ogive-Cylinders. 


Figure  13.  Calculated  Values  of  Nose  Pressure  Drag  Presented  as  a 
Function  of  Nose  Length  and  Nose  Bluntness  for  Varying  Mach  Numbers 


Nose  Pressure 


rapidly  increases  with  decreasing  nose  lengths  when  these  lengths 
are  less  than  1.5  calibers.  As  the  body  becomes  more  blunt,  the 
flow  expands  more  rapidly  and  a larger  pocket  of  supercritical  flow 
is  formed.  It  is  this  phenomena  which  causes  the  nose  drag  increase 
with  decreasing  nose  length. 

The  nose  drag  coefficient  versus  nose  radius  with  nose  length 
as  a parameter  is  presented  in  figures  13b,  d,  f,  h.  These  results 
for  each  Mach  number  are  obtained  by  plotting  the  results  of  the 
faired  curves  in  figures  13a,  c,  e,  g for  a constant  length.  This 
method  enables  the  smoothest  set  of  curves  from  the  least  amount  of 
calculated  data  to  be  obtained.  The  final  coefficient  values  from 
these  figures  are  presented  in  Table  I.  This  table  is  then  incor- 
porated into  subroutine  TRANS  of  the  NSWC  computer  code. 

Figure  14  gives  the  wave  drag  obtained  by  solving  the  transonic 
perturbation  equation  and  the  wave  drag  obtained  by  solving  the  full 
Euler  equations  for  various  tangent  ogive  nose  lengths  throughout 
the  transonic  Mach  number  range.  The  new  results  extend  the  lower 
end  of  the  nose  length  curves  from  1.5  calibers  to  0.75  calibers. 

To  test  the  accuracy  of  the  new  method  for  calculating  the 
nose  drag,  a suitable  set  of  test  cases  were  run  and  compared  with 

9 

an  available  experiment.  These  cases  involved  variations  in  nose 
bluntness,  nose  length,  and  transonic  Mach  numbers  and  are  presented 
in  figures  15-17. 

In  all  three  test  configurations  the  afterbody  is  10  calibers 
in  length  and  the  nose  bluntnesses  were  0.0,  0.250,  .50, 

and  .750  calibers. 

The  following  observations  can  be  made: 

1.  For  the  pointed  tangent  ogive  body  the  revised  computer 

program  resulcs  compare  more  favorably  with  experiment  than 
the  previous  computer  program  results.  This  improvement 
is  as  high  as  39  percent  over  the  previous  result  for  a 
pointed  2 caliber  nose  length  projectile  at  a free  stream 

M„  = 0.8. 
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a)  Two  caliber  tangent  ogive  body 
with  a pointed  nose. 
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b)  Two  caliber  tangent  ogive  body 
with  a nose  bluntness  of  .125. 

Figure  15.  Total  and  Forebody  Drag  Comparisons  of  the  Present 
and  Previous  Theoretical  Results  with  Experiment 
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c)  Two  caliber  tangent  ogive  body  with 
a nose  bluntness  of  0.25 


d)  Xwo  caliber  tangent  ogive  body  with 
a nose  bluntncss  of  0.375 


Figure  15.  Concluded 

i i 


36 


a)  Three  caliber  tangent  ogive  body 
with  a pointed  nose 
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b)  Three  caliber  tangent  ogive  body 
with  a nose  bluntness  of  0.125 


Figure  16.  Total  and  Forebody  Drag  Comparisons  of  the  Present  and 
Previous  Theoretical  Results  with  Experiment 
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a)  Four  caliber  tangent  ogive  body 
with  a pointed  nose 
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b)  Four  caliber  tangent  ogive  body 
with  a nose  bluntness  of  0.125 
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Figure  17.  Total  and  Forebody  Drag  Comparisons  of  the  Present  and 
Previous  Theoretical  Results  with  Experiment 
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2. 


For  the  larger  nose  bluntnesses  the  revised  computer 
program  results  again  compare  more  favorably  with 
experiment  than  the  previous  results.  The  maximum 
improvement  of  27  percent  occurs  at  a free  stream 
= 0.8,  a nose  bluntness,  = .375,  and  a nose 

length  L^/2R^  = 4. 

3.  For  transonic  Mach  numbers  less  than  0.85  the  revised 
computer  program  gives  results  which  are  closer  to 
experiment. 

4.  The  experimental  points  were  obtained  in  laminar  flow 
whereas  the  calculations  are  performed  assuming  turbu- 
lent flow.  This  could  explain  the  discrepancy  between 
the  theory  and  the  experiment. 
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There  have  been  no  changes  in  the  input/output  portions  of 
the  NSWC  code.  The  changes  in  TRANS  and  MAIN,  and  the  additional 
subroutines  necessary  are  documented  below. 


INTERD,  INTERS 

These  subroutines  are  used  to  interpolate  for  the  value  of  a 
function  at  a given  point.  Three  point  Lagrange  interpolation  is 
used . 
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The  changes  in  TRANS  are  shown  by  a solid  line  down  the  left 
hand  margin.  A sample  input  data  sheet  is  shown  in  figure  18 
where  an  accompanying  sketch  with  the  geometric  parameters  describing 
the  configuration  is  also  given.  Figure  19  gives  the  sample  output 
for  the  test  case.  for  0.5  < < 0.8  is  reduced  quadratically 

to  zero  from  a table  value  at  = 0.8.  This  and  other  improved 
modifications  to  program  changes  in  the  original  NEAR  TR  153  were 
made  by  L.  Devan  at  NSWC. 

CONCLUSIONS  AND  RECOMMENDATIONS 

1.  An  accurate  implicit  axisymmetric  Euler  equation  solver 
has  been  used  to  obtain  nose  drag  coefficients  for  transonic  flow 
over  sphere  ogive  cylinders. 

2.  A matrix  of  drag  values  was  calculated  for  the  following 
parameters 

0.8  < M„  < 1.2 

0 i R^/D  i . 5 

.75  L /D  10 
n 

3.  The  computed  matrix  was  incorporated  into  the  NSWC  computer 
code  as  a table  lookup.  Several  test  cases  were  run  and  compared 
with  experiment. 

4 . It  is  recommended  the  Euler  equation  solver  be  used  to 
calculate  drag  coefficients  for  bodies  with  surface  discontinuities, 
and  be  used  to  calculate  drag  coefficients  for  sphere-ogive  bodies 
in  the  Mach  number  range  0.5  to  0.8. 
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0ATA(AM(T  1,1  = 1.41  /. 9, .95, 1.05, 1.2/ 

OAT  A ((CA15( I . J) , T = 1 .41  . J=1 , 10)/.  4 22. . 51 8 , . 6 7 , . 82 , . 2“  , . 4 , . 541 , . 652, 

1 .198. . 311. .441. .552..  130..250,.  3 65.  . 497, .0  89,. 211,. 3 35,. 47,. 065, 

1 .19.  . 31. . 451 . .06. .18. .291. .435,  . 058, .18, .285, .425  ..06, .184,. 3, .435 
),. 38, .278  , .41, .565/ 

RATA ((CA2  (I ,J) ,1  = 1,4)  ,3=1,10)/.  29,.  395 , . 5 0 5, . 64, . 2 , . 295 , .409, . 521 

1. . 13.. 22 3.. 3 3 9.. 4 5,. 0 81, .175,. 29, . 43 4,  . 051 ,. 144, .26 ,. 374, .04, . 13, 

1 .24,. 3 5, .0  4. .123, .229, .335, .045, . 128 , .22  3, . 333, . 0 55 , .15, . 24, .355, 


1 
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^•b6t  *2799 f 9 / 

)f«0ASfal5««37t«37**C*>lt9ll5«*23« •333t«93t*091f«202t*302«*0?«*PAt 
l•lAf•291«•022«•&81t•X7««27t•C31«•091«•12«•2A««Q42f•121f•2t•32«•0A« 

1 .278. . 41. .5d9/ 

OATfl  t «C4  3 (I .J»,I  = 1.4> .Jr  1, 107.151,  .198,.  3 02, .353. . 101,. 15..25, 

1.319.. 3b5..11..21,.29,.C3-,.C  75  ..18,  .2S3,.02,.0  55 ,.16, .245, .018, 
l.u47,. 145, .231, .022, .049,. 1 35, .23  ,.3 31.. 065  ,.13, .24,  .042,. 13 8, 
>.175, .30 2, .08, .278, ,41, ,5657 

nA7A<(CA4  (I,J),I=l,4»,Jrl,lC)7.1,,130,  . 2 1 1 . .225 , . 07 0 , . 091, . 1 75 , 
I.20  5..04,.361,.14,.l9,.a2..0  42,  .115,.17;,.0i,.0<,,.i,.17,.01,  .031, 

). 095, . 17. . 02, . u 48, .1, .16, .0  31, .06 5. .121, .211, .042. .10*, .175, .3  9, 
>.O8,.£70,.*,1  ,.5657 

OAT  A t (CA6  ( t , J)  , 1 = 1,41  , Jr  1,  lOt  7 .068,  .080,  .105,  .095,  .041,  .049,  .085, 

> .092, .02, .038, .067, .095, .007, .032, .055, .1, .006, .031, .051, .109, .01, 

I .031, .061, .121, .02, .048, .082, .151, .031, .065 , .12, . 205, .042, .108, 
).l75f.29,.ij«,.27«,.41  ,.565/ 

04T4I(C»8  (I.J>tI=lt4),Jrl,lCI7.35,.04,.34,.04,.029,.031,.335,.04, 
1 .012, .021, .03, .05, .007, .028, .03, .06, .006, .028, .038, .07 , .01, .031, .05, 
).ia5«.U2t.S4e..v7X.145,.  3*1, .065, . 12  > .23  4t .04 2, .10*, .175,. 29, .08, 

1 .278. . 41 . .5557 

nAT4Mf;ii3(i,j»,r=i.4>,jri,ij)7  .041,. 025,, 025, .025,. 02, .020, .025,. 025, 

) .008, .015, .025, .03, .007, .020, .025, .041, .006, .021, .025, .055, .01,, 031 , 

>.G4, .095, .02,. 049,. 07,. 14, .031, . 363, .12, .204, .042, .108,. 175, ,?9, 
>,u8,,278,.41 ,.5657 
C8D*0. 

IF199P1NN).LT.-C.uajll  GO  TO  "7 
GO  TO  98 

87  V0\»sV0V5-.S498 

IFCVOV.LT.l.  ) GO  TO  89 
ARe'^  = 3.14159*RRE*^»*2 
tF{RSH4i’6.Lr.4)  GO  TO  88 
16(9911.60.2)  GO  70  1 
I61N9L')NT.E0.2)  go  to  1 
JrN92tl 
GO  TO  2 
1 Jr'I'IS*! 

2 00  13  LrJ,N9 
XXr<9(L) -X0( J) 

OELTArlTAMd./  t^.•l^lLt  I 

I6(R9P  tJ-31.l.T.R3P(l>>  06LT  ftr  ftTA  9 ( . 2 7 AXL) 

CA9Ari.4 
Cl=l .♦Gl«l 
00=5007(01) 

C3=OOV  ••2 

C4ri .-C3 

05=047 (C1*C3 I 

C6  = 3 .*051747 (2. *001 

07=25. •C1*V0V  **(2.73.) 

C8r.5*C47(Cl*C3l 

C9=1.25*C5**2 

C13=2.*C57(00W  **(2.73.) )*C6**(2.73. ) 

Cll=(C67V0V  )**(4.73.) 

CS1rC7*{CS*(C9*Clj*011)**(.5) ) 

C*S09T (C50) 

T=2.*A14 

CPI  =.4* (r-C)750RT (C1*V0V  •* (2.  73. ) I • ( .04*(V-r)**27CCl*WOV  •* 

1(2.73.)) -C47(C1*C3))**.5 

I6(t.gt.c)  :pi=t. 

0ELTAr-RBO(L ) 

0494=1.4 

C1=1.*G494 

Cfa»SOPT(Cll 
C3='(0V  **2 
C4*l .-C3 
C5»C47 (C1*C3  ) 
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ooo  oooooo 
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JBOM  OQPY  EURBlSHIiD  TO  HOC 


Cb«i  .‘OtL  I A/  Jl 

C7»25.*C1*V3 V **(2.73.» 

C8*.S*C<./  (Cl  *03) 

Cq=1.25*C5**2 

C10*2.*C5/(VOV  •*(2./3. I l*C6»*(2.  /J.  I 
C11*(C5/V0V  )••(<*. /3. I 
CS0=Cr*(C8»( C9*C10*C1H  *♦(.51) 

C=S3RT(CSQ) 

r=)(X*2. 

CPWIL,  ll=.l»*IY-C»/S1»T  (Cl*VOV  • • ( 2.  /3  . I ) • ( . 0 <••  ( V-C  ) **2/ (Cl*  *0V  •• 
l(2./3.)» -C4/ (C1*C3»)** .5-nELTA**2*:»l 
IF(Y.GT.C)  CPV(L,1I=CP1 
15  00  10  Kl>1.7 

CPY(L,Kl)=C»V(Ltl) 

10  CONTINUE 
JA<J 
JB«NN 
SUN1*0. 

SUN2=0. 

$UN3=0. 

CALL  SIMP 
CB3«2. •SUMl/AREF 
VO-YOV-1  . 

IFIYO.LE. 0.061  CBlsCMO 
V1*V0V$ 

89  CONTINUE 
IFfYOYS.GT.O .95»  GO  TO  90 
C93=0. 

GO  TO  88 

90  ceO^CPl* {V0YS-.95)/  (V1-.95I 
88  CONTINUE 

C TABLE  LOOKU®  IS  9ASF9  ON  LN/RR,  RN/?9 

ANL1>ANL 

ANLsANL*RR*(-1.»1./(1. ♦SIN(THU>  » 

ANL=2.*ANL 

VVaYOVS 

IFfVOVS.  LT.  .8»VOVS=.8 


DETERMINE  RNOR8,  THE  CORRECT  RA0IU3  OP  THE  NOSE,  TO 
BE  USED  IN  TABLE  LOOKUP. 

RN3R9=2.0*R( 1) /COS(THl) 


HRITE(6,2aOI  VOVS,ANL .RNORB 

200  FORMATC*  *,*MACH  NO.,  NOSE  LENGTH,  RN/ 29 » , 3 F15 . 10, /// 
1 ’AS  ANO  DRAG*) 

IFIANL.GT.IC.)  ANL=ia. 

IF<ANL.  LT.  1.5)  ANL=1.5 
IFfRNORB.GE.fl.M)  GO  TO  21 
Jl»l 
J2  = 5 

GO  TO  22 

21  IFIRNORB.LE. 0.51  GO  TO  23 
Jl«5 

J2310 
GO  TO  22 
23  Jl»i* 

J2«8 

22  DO  16  J«J1,J2 
no  13  1x1,6 
IF(ANL.Gr.6.  ) GO  TO  26 
CA15R(T)*CA15(I,J) 


UAi  •<ii»=CRi; 

CA25R(r)=CA’5(I,J) 

2«*  CA3‘?{II  = CA3(  r.  Jl 
CA4  R(T>=CA4  (I,J) 

IF(ANL.LE.4. ) GO  TO  1 3 
CA6  R(I)=CA6  »I,J> 

CAS  R(II=CAS  (I.J) 

CAljR(II=CAlOCI,J) 

13  ''ORTtRUE 
IFUNL.GT.4.)  go  to  2f) 

19  CALL  INTeR5{RN0R9,  .3.  .4,  .9,.6,.r,C)(l»).Ca(SI,CD(6l,Cn(n,C3|St, 

) CAN) 

CALL  INTERD( AM,CA19K, tfOVS,AO< J) t 4t2» 

CALL  INTEROC AM,CA2  R. VO VS , A1 < J) . 4 , ») 

CALL  IMTEPOC  AM,CA25R, VOVS, A2 « J» t 4, >) 

26  CALL  INTERO( AM,CA3  R. VO VS , A3 « J)  , 4 , ») 
call  INTEROr AM,CA4  R, VOV S . 44 f J) , 4, ?) 
tF(ANL.LE.4. » GO  TO  26 
CALL  INTER0»AM,CA6  R t VO VS . A5 ( J) , 4 , 2 » 

CALL  INTERO* AM.CAS  R, VOV S , A6 ( J) , 4 , » » 

CALL  INTEPD<  AM , C A lOR, VO VS , A 7 ( J) , 4 , 2) 

GO  TO  12 

25  call  INTER5( ANL,1.5,2. .2.5,3.,4.  . A3IJ).A1(J),A2(J>, A3(J),A4( JI, 

I CO(JI) 

GO  TO  14 

12  CALL  INTER5( ANL.3..4. ,6. ,S. ,10.  . A3( J) .S4( J) , A5 (J) ,SeiJ) ,A7(J), 

) CO(J)» 

14  CONTINUE 

IF(RNORO  ,GE.  0.4)GO  TO  16 

CALL  INTER5(RNORO,0.,.1,.2,.3,.  4,CH1),C0C2)  , CO  ( 3 ) , COI 4)  , COt  5»  , 

) CAN) 

GO  TO  17 

16  IFCRNORR  .LE.  0.5)GO  to  19 

CALL  INTFR5C  RNORR,  .5,  . 6 . . 7,  . S ,1 . , C ) ( 6) . C 0 ( 7 ) , CO  ( 8 ) , CO  ( 9)  , CO  1 10  > , 
I CAN) 

GO  TO  17 

17  IF(VV  .GE.. 8)GO  TO  IS 
CAN=CSN* ( <VV-.5)/.3) **2 
vovs=vv 

18  CAN=CAN+COO»CAH 
NPITE(6,215)  CAN, CAW 

215  FORMATt*  •,»NOSE  DRAG  ANO  HAVE  ORAJ* , 2ti 5. 1 0 I 
ANL=ANL1 
99  return 
ENO 

CHANGE  OCCURS  AT  HAVE. 8 

-TeMMON^OrPRyXir  tRH.OiOEF — 

C0MM0N/GE04/K,F,RR,RRFP,thi 


IBIS  PA®  IS  BfiST  QUALITY  PRACTlCABLt 
KtOIIOOlPY  FURNISHED  10  QDC  — ^ 


47 


j|i 

■■■■■■■■■■■■■■■■■■■■■■■■I 
■■■■■■■■■■■■■■■■■■■■■■■■I 

K^-aairiaaaaBaaaai 

irtBa?knaaaaBaai  

IClMKiaBHHBBHHBBBBnBagi 
IBBBL*  BBBBBBBBBBBBgBBaaBl 
ICTBBBBBBBBBBBBBr 
IkiBBBBBBBBBBBBBI 
IL*  ■■■■■■■■■■BBBI 
lO^BBBBBBBBgBBBBI 
ItSiBBBBBBBBiBBBBI 
KsiBBiBBBBBaaggg! 

!igsr:s:::B:s:&B.BBHB..Bi 

l■BBB■BBBa■BBgBB■■■■■■■i■■! 
■aaai*;aaaaaaaaaaar 
HSBBGTBBBBBBBBBBBI 
Koaaaaaaaaaaaaaai 
lljSBBBBBBBBBBBaaai 
■s<  aaaaaaBBBBBBBBi 
■s  aaaaaaaaBBBaBBi 
■Maaaaaaaaaaaaaai 
K^BBBBBaaBBaBBBBI 
■■aarraaaaaaaaaaai 

■SaaSgBaaaaBaaaai 
■■BBgBBBBBBBBBBBI 
■aBBBBBBBBBBaBBBI 
■BBBaBBaBBaBBBagi 
■■I2BSBBBBBBBBBBBI  

ISBBala^aBBBBaBaBBBaSaBaSSI 

ll^BBBBBBBBBBaBaBBBBBSBBBBBI 
■■BBL^BBBBBBBBBBBBBBBBBBBBgi 
■CiBBBBBBBBaBBBaaaBBBBaaaBBI 

lEB^BBBaBBgaBBBBBBBBBBBai 

SggSBSSBBBiSSS&BSBaSBBaSS 

■BEi  • bBBBS 

■BaBKSSBBaBBBBBBBBBBBBBBBBi 
■BBBBBBBBBBBBBBBBBBBBBBBBBL. 
■BBgBBBBBBBBBBBBBBBBBBBBBBBl 
IBBUBBBBBBBBBBBBflBBBBBBBBaal 
■L^UKaaBaBBaBaBaBBagaBBagagpi 
IBBBtnBBdBBllLBQBBBLIiBBBBBBBiiBrLl 
■BBB^BB:'i;LLSiii:ssL^:«:BS!-:;MB^pL^^[2^r:« 

■BBBBI 

■Bl^iSB' 

■BBB^i 

■aai'L^i 

■BBL*JBBBBBBBBBBBBI 
■BBBaBBBaBaBBaaai 
■CI!.CBBBBBBBBBBBBI 
■BBBli:  BBBBaBBaBBaBBBBBBBBBBI 
■aaaaBBBBBBBBBBBBBBBBBBBBBI 
■BBGUBBBBBBBBBBBBBBBaBBBBBgi 
■BBBBaBaBBBBBBBBBBBBBBaaBBI 
■BL.*MBaBaBBBBgagBBBBBBBBBBBI 

■BBB:7'BB:?vii:B.';rii»m'ii£a(ciii'izBC?  i'iizl. 

■ BBBaBBBKBBKBBBBBBBKiaBBBEBBl 

IL»ll^CBB>dL»;BBBBBaeCSgeceC£g.{Iggl 


CA»t  NU,  1 

,or  ATTACK  « .04000669 

reference  diameter  « I.OOOFT 

R6FER6NC6 

CUNDITIUnS 

SPIED  OF 

SUUNU 

« 1116,890  FT/StC 

- - PENSITV 

■ ,0023769  SLUCS/FTJ 

ABSOLUTE 

VISCOSITY 

■ ,000000374528  L8»StC/FT2 

MOMENT  REFERENCE  • 

O.OOCALIBERS  FWU^‘  NU8E  TIP 

- 

BOOT  COORDINATES 



R 

0,0000 

0.0000 

- - .U50 

.0050 

• 2500 

.0800 

. 

.3750 

.1200 

4?>000 

.1550 

• 6250 

.1900 

• 7500 

.2250 

_ ,8750 

.2500 

1,0000 

.2800 

1.1250 

.3100 

- . 1.2500 

.3350 

1,37S0 

.3600 

1,5000 

.3800 

1.6250 

• 4000 

1,7500 

.4200 

1,8750 

.4500 

_ . 2,000.0 

• 4500 

2,1250 

.4600 

2.2500 

.4750 

2.3750 

.4800 

• 

2.5000 

.4650 

2.6250 

.4900 

- _ 2,7500 

,4950 

2.8750 

.4975 

3,0000 

.5000 

- -4.0000,. 

• 5000 

6.0000 

• 5000 

- ■ . . .. 

8.0000 

.5000 

10,0000 

• 5000 

. 

Figure  19.  Sample  Output  for  Tangent  Ogive  Configuration 

of  Figure  18 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLa 
ttKm  COPY  rURMlSHED  TO  WC  


BODY  AXUL  ^URCL  CONTKlHUT  lUNS 


hACH  nO.OkI*^ 

XHICTIUS  base 

PRESSURE 

PRESSURE  ■protrusions' 

• 

total 

1.500 

.Uftii 

,1910 

• 1504 

0,0000 

,4046 

l.iOO 

,0601 

,2100 

,1524 

0,0000 

• 43S4 

1,050 

,0700 

,2190 

• 1151 

0,0000 

,4048 

1,000 

,0717 

.2090 

• 096i 

0,0000 

,3770 

,050 

,0726 

.1044 

,0660 

u,ouoo 

,3231 

,eoo‘ 

.0755 

• 1350 

j0767 

;,oooo 

*2e72 

,000 

• 07»§ 

• 1205 

.0085 

0,0000 

.2168 

total  static  AtRO0YKAMlCS(f QHCt/ALPHA) 


ACH  NO,  CO 

CN  Cl 

CN 

CNAl 

CRAl  XCP/u 

1,500  ,4046 

0,0000  0,0000 

0,000 

0,000 

0,0000,0000 

1,200  ,4184 

O.OOOOO.QOOO 

0,000 

0,000 

0,0000,0000 

1,050  ,4048 

0,00004,0000 

0,000 

0,000 

0,0000,0000 

1,000  ,3770 

0^00000,0000 

0,000 

0,000 

Q.OOOQ.QOOO 

,950 .1231 

0,00080,0000 

0,000 

0,000 

0,0000,0000 

,800  ,2872 

0,00000,0000 

u,000 

0,000 

0,0000,0000 

,600  .2083 

O.OOOOO.OOCLO 

0,000 

0,000 

0*QOO0«000O 

MACM  NO, 

CLP 

MWiir  # 1 A 66'^  V 1 '.vrt  El  A W ^ .f  cr.i  4 r.tmW 

CM06CRA0 

CNPAl 

1,50 

•,042 

•369,240 

0,000 

1,20 

•,042 

•369,240 

o.»QOo  ; 

1,05 

••04} 

•369,240 

0,000 

1,00 

• ,044 

•369,240 

0,000 

• 95 

• .046 

•169^240 

OaOOO 

tto 

• .066 

•369,240 

o.ooe 

.60 

••066 

•369,240 

0.000 

CNHA5 


0.000 

o.ctoo 

0,000 

0.000 

Q«000.' 

o.soo 

0.000 


Figure  19.  Concluded 
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GLOSSARY  OF  TERMS 


A Jacobian  matrix  defined  in  Equation 

B Jacobian  matrix  defined  in  Equation 

1/2 

c speed  of  sound  = (Yp/p) 

Cjj  drag  coefficient 

Cp  pressure  coefficient 

dy  space  increment 

dt  time  increment 

D drag 

e energy  per  unit  volume 

E matrix  of  conservative  variables  defined  in  Equation  1 

E + CyF)/J 

F matrix  of  conservative  variables  defined  in  Equation  1 

F (n^u  + + nyF)/j 

H source  matrix  defined  in  Equation  1 

I identity  matrix 

j j = 0 2-D  flow;  j * 1 axisymmetric  flow 

J Jacobian 

length  of  nose 

M free  streeun  Mach  number 

p pressure 

1 2 

q dynaunic  pressure  = 

2 

q total  velocity  vector  squared 

rb»Rjj  radius  of  the  cylindrical  afterbody 

r ,R  radius  of  the  spherical-  nose 
n n 

S smoothing  term  defined  by  Equation  7 

2 

S^^j  reference  area  * irr^j 


A 3E 

4.  = — 

3u 
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t time 

u velocity  component  in  x-direction 

U matrix  defined  in  Equation  1 

u U/J 

V velocity  component  in  y-direction 

X coordinate  direction  in  physical  plane 

y coordinate  direction  in  physical  plane 

yj^  radius  of  the  body 

Y term  included  in  source  term,  H of  Equation  1 

y ratio  of  specific  heats,  1.4  in  this  study 

6 central  difference  operator  in  Equation  6;  shock  angle  in 

figure  8 

e smoothing  coefficient,  0 i e ^ 0.4 

n coordinate  direction  in  computational  plane 

C coordinate  direction  in  computational  plane 

p density 

a eigenvalue  defined  in  Equation  5 

T time 

0 meridional  angle  defined  in  figure  2 

Subscript 

b body 

1 increment 

j ^-direction 

k n-direction 

n,N  nose 

Superscript 
n time  direction 

( ) averaging 
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